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Abstract. 

We discuss aspects of open and hidden charm production in hadron-nucleus 
collisions at RHIC and LHC energies. We hrst discuss the extraction of the total charm 
cross section in lower energy collisions and how it compares to next-to-leading order 
quantum chromodynamics calculations. We then describe calculations of the transverse 
momentum distributions and their agreement with the shape of the measured STAR 
transverse momentum distributions. We next explain how shadowing and moderate 
^C) ■ nuclear absorption can explain the PHENIX J/ifj dAu/pp ratios. 
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detected, completing the picture of the fourth quark begun when the J/ip was detected 
in pBe and e + e~ interactions. The charm quark was postulated to have a mass between 
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1. Open charm production at RHIC 

Open charm measurements date back to the late 1970s when D and D mesons were first 



1.2 and 1.8 GeV, within the regime of perturbative quantum chromodynamics (pQCD). 
Because of its relatively large mass, it is possible to calculate a total cc cross section, 
not the case for lighter flavors such as strangeness. Charm hadrons are usually detected 
two ways. The reconstruction of decays to charged hadrons such as D° — > K~n + (3.8%) 
and D + — > K~-k + 'k + (9.1%) gives the full momentum of the initial D meson, yielding 
the best direct measurement. Charm can also be detected indirectly via semi-leptonic 
decays to leptons such as D — > Klv\ although the momentum of the parent D meson 
remains unknown. Early measurements of prompt leptons in beam dump experiments 
assumed that the density of the dump was high enough to absorb semi-leptonic decays 
of non-charm hadrons, leaving only the charm component. At modern colliders, it is 
not possible to use beam dumps to measure charm from leptons but, at sufficiently high 
Pt, electrons from charm emerge from hadronic cocktails pflE]- 

Although D mesons are usually used to determine the total cc cross section, other 
charm hadrons also exist. The excited D states, D*s, decay primarily to charged and 

* This work was supported in part by the Director, Office of Energy Research, Division of Nuclear 
Physics of the Office of High Energy and Nuclear Physics of the U. S. Department of Energy under 
Contract Number DE-AC03-76SF00098. 



2 



neutral D mesons. The charm-strange meson, the D S} decays to charged hadrons and 
simultaneously to leptons. The lowest lying charm baryon is the A+ which decays 
primarily to A(uds) but also decays to pK~ir + and to the lepton channel. The heavier 
ground state charm baryons and their excited states (E c and higher) decay through the 
A c channel. The charm-strange baryons are assumed to be a negligible contribution to 
the total cross section. A selection of charm hadrons, their masses, decay lengths and 
branching ratios to leptons and charged hadrons are given in Table ^ 
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Table 1. Ground state charm hadrons with their masses, decay lengths (when given) 
and branching ratios to leptons (when applicable) and some prominent charged hadron 
decays. 



Extracting the total charm cross section is a non-trivial task. To go from a 
finite number of measured D mesons in a particular decay channel to the total cc 
cross section one must: divide by the branching ratio; correct for the luminosity, 
od = No/ Ct; extrapolate to full phase space from the finite detector acceptance; divide 
by two to get the pair cross section from the single Ds; and multiply by a correction 
factor [Hj to account for the unmeasured charm hadrons. There are assumptions all 
along the way. The most important is the extrapolation to full phase space. Before 
QCD calculations were available, the data were extrapolated assuming a power law 
for the xp distribution, related to the longitudinal momentum of the charm hadron by 
xp = p z /(VS /2) = 2rriT smb. y/\/S. The canonical parameterization is (1 — xp) c where c 
was either fit to data over a finite xp range or simply assumed. These parameterizations 
could lead to large overestimates of the total cross section when < c < 2 was assumed, 
especially when data were taken only near xp = 0. Lepton measurements were more 
conservative but were typically at more forward xp. 
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1.1. Total cc cross section 

There has been a great deal of improvement over the last 10-15 years. Next-to- 
leading order (NLO) calculations are used in the phase space extrapolation, resulting in 
considerably less ambiguity in the shape of the x F distributions, da/dxp. The transverse 
momentum distributions are more difficult, as we will discuss later. To calculate the 
total cross section to NLO, scaling functions [3] proportional to logs of p 2 /m 2 are useful 
where /i is the scale of the hard process. The hadronic cross section in pp collisions can 
be written as 

a pp (S,m 2 ) = J dxidx 2 fi(x 1 ,/j, F ) fJ(x 2 ,i^ F )a ij (s,m 2 ,ij, F ,iJ 2 B ) (1) 

where x± and x 2 are the fractional momenta carried by the colliding partons and ff are 
the proton parton densities. The partonic cross sections are 

m 2 
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with s the squared partonic center of mass energy, p = 4m 2 /s and f^' 1 ^ are the scaling 
functions given to NLO in Ref. It is most consistent to assume that the factorization 
scale, (xp, and the renormalization scale, /ir, are equal, fi — \ip = /ir. There is no 
dependence on the kinematic variables. Some NNLO calculations are available near 
threshold, s = x^ 2 S ~ 1.3 (4m 2 ), applicable only for y/S < 20 - 25 GeV H U 
The NLO corrections to the leading order (LO) cross sections are relatively large, 
K^ x > = 0"nlo/°"lo ~ 2 — 3, depending on /i, m and the parton densities [S]. The NNLO 
corrections are about as large to next-to-next-to-leading logarithm |lj but decrease to 
less than K^ l > when subleading logs are included j^j- This K factor is large because, in 
the range 1.2 < m < 1.8 GeV, m < fi < 2m with a 5-flavor QCD scale, A 5 , of 0.153 GeV 
for the GRV98 HO and 0.22 GeV for the MRST parton densities, 0.21 < a s (c) < 0.4, 
nearly a factor of two variation. (The larger value corresponds to the smallest m and 
fi values with the larger A5.) The situation improves for bottom where a s is smaller, 
0.16 < a s (b) < 0.28, and is quite good for top, 0.092 < a s (t) < 0.12. Instead of 
presenting a wide range of possible cross sections and emphasizing the uncertainties, 
the approach taken in Ref. [Z] has been to "fit" m and \i for a particular parton density 
and extrapolate to higher energies. The results are compared to some of the total cross 
section data [8^ on the left-hand side of Fig. ^ The data tend to favor lower values of 
m, 1.2 — 1.3 GeV. The two curves cross each other because the MRST calculation with 
/1 = 2m increases faster at large \fS and smaller x due to the stronger QCD evolution 
of the parton densities at the higher scale. Although the fixed target results are in 
good agreement with the calculations, the PHENIX point [9J at 130 GeV, from Au+Au 
electron measurements, and the STAR point ^Ul, from a combination of electron and 
reconstructed D measurements, are generally above the calculations. The STAR point 
is about a factor of four over the calculated cross section. The higher energy pp data 
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from UAl jTTj and CDF ^2] are in better agreement with the calculations. (At these 
energies, the pp and pp cross sections differ by less than 1% for y/S > 630 GeV.) 
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Figure 1. The NLO total cc cross sections as a function of \fS (left-hand side) and 
charm quark px distribution at y/S = 200 GeV in the range \y\ < 1 (right-hand side) 
in pp interactions. The curves are MRST HO (solid) with m = 1.2 GeV and fi 2 = Am 2 
and GRV98 HO (dashed) with m = 1.3 GeV and [i 2 = m 2 . 



1.2. Open charm transverse momentum distributions 

Now we turn to the transverse momentum, pr, distributions. In this case, the quark 
mass is no longer the only scale and pr-dependent logs also appear. Thus, to interpolate 
between the high p T scale of px and the low px scale of m, a scale proportional to m T , 
the transverse mass, is the natural choice. The charm quark p? distributions are not 
strongly dependent on quark mass for p? > 3 GeV, as may be expected, where the 
difference in rate is ~ 20% between m = 1.2 and 1.8 GeV. The difference in the total 
cross sections is almost all at px < 3 GeV. Changing the scale changes the slope of the 
Pt distributions. The distributions are harder for /x = m than \x = 2m. The average p?, 
(pr), increases with m and is larger for \i — m. 

More modeling is involved for D mesons in the treatment of fragmenta- 
tion/hadronization and momentum broadening. If factorization holds in the final state 
(universal fragmentation functions) as well as in the initial state (universal parton dis- 
tributions) then the fragmentation functions extracted in e + e~ should also be applicable 
to pp and dA. However, this assumption does not work well for charm. The Peterson 
function, generally used in hadroproduction codes, reduces the charm hadron momen- 
tum by 30% relative to the charm quark. As shown in Huang's talk ^H], the Peterson 
function agrees reasonably well with the e + e~ data. (However, it does not include any 
scale evolution. In low v^S* collisions, the momentum reduction due to fragmentation 
can be compensated by intrinsic transverse momentum, fcy, broadening. However, such 
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broadening cannot compensate the xp distributions, only marginally affected by kp 
smearing. We have previously shown that the D meson xp distributions are consistent 
with no momentum loss during charm quark hadronization [14].) The exclusive NLO 
QQ code of Ref. [T3] includes fragmentation and broadening. This program adds the kp 
kick in the final, rather than the initial state. The initial kp of the partons could have 
alternatively been given to the entire final-state system, as is essentially done if applied 
in the initial state, instead of to the QQ pair. The Gaussian function g p (kp), 

9 P ( k T) = } 2 \ exp(-4/ (kp) P ) (3) 

7l\Kp)p 

multiplies the parton distribution functions, assuming the x and kp dependencies 
completely factorize. If true, it does not matter whether the kp dependence appears in 
the initial or final state. There is no difference if the calculation is LO but at NLO an 
additional light parton appears in the final state. The difference in the two methods is 
rather small \ik\<1 GeV 2 [8J. The value (k\) p = 1 GeV 2 was used in Ref. 0. 

The effects of fragmentation and intrinsic kp broadening of (k^,) = 1 GeV 2 
compensate each other at \/S = 20 GeV to give a D meson p T distribution very similar 
to that of the charm quark [T^j. However, at RHIC energies, the situation is quite 
different. Due to the higher (jpp) at larger \/S, the effect of broadening is relatively small 
and cannot compensate for the momentum loss induced by fragmentation. Interestingly 
enough, the STAR D and D* pp distribution agrees rather well with the calculated NLO 
charm quark distribution, as shown in van Leeuwen's talk ^Ej. On the right-hand side 
of Fig. [TJ we show the corresponding p T distributions at \/S = 200 GeV for the two sets 
of parameters in the total cross section curves on the left-hand side. The differences 
in the slopes are due to the different scales while the normalization difference is due 
to the choice of charm mass and the parton densities — the MRST densities generally 
give a larger cross section due to their larger a s . However, the curves need to be scaled 
up by a factor of four to agree with the STAR normalization ^H], as may be expected 
from the total cross section results. The shape of the charm quark pp distribution at 
y/S = 1.96 TeV also agrees quite well with the CDF data from the Tevatron ^2]- Given 
the large discrepancy between the pQCD result and the STAR cross section, it might be 
surprising that the normalization is also in good agreement with the sum of the charged 
and neutral D data scaled to include D s and A c production. No total cross section is 
available because only charm hadrons with pp > 5 GeV have been measured so far. 

Other model calculations of charm production at collider energies are available. 
The FONLL calculation resums logs at pp 3> m, resulting in a harder charm quark 
distribution and uses a correspondingly harder fragmentation function to get agreement 
with the CDF data. It underestimates the low pp STAR data, resulting in an even 
lower total cross section than NLO pQCD JH]- A calculation with unintegrated gluon 
distributions and /c^-dependent matrix elements, assuming saturation behavior at low 
x, has also been made [HIJ- However, the x values of the STAR data are not really very 
low. At RHIC, from kinematics x ~ 0.01 at y = 1 and pp = 0, the highest rapidity 
measured by STAR by kinematics alone. In reality, the actual (x) may be higher when 
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weighted by the parton densities. At higher p T , x is larger still, suggesting that the 
applicability of small x physics for charm at RHIC is rather dubious. 

Finally, we would like to discuss reasons why fragmentation does not seem to 
factorize for charm, as expected. Factorization breaking has been suggested from studies 
of the xf distributions of e.g. D + and D~ production, particularly in tt~A interactions 
where the D~ is leading relative to the D + since the D~ shares a valence quark with the 
7r~ while the D + does not. Several mechanisms such as intrinsic charm [2] and string 
drag have been proposed, both of which involve charm quark coalescence with spectators. 
Such comoving partons are naturally produced in a hard scattering. Although it is not 
intuitive to expect coalescence to work at high p?, it seems to do so for charm. 

2. Nuclear dependence of J /if) production at RHIC and LHC 

We now turn to J /if) production in d+Au interactions at RHIC. Previously, we calculated 
the effect of shadowing alone on the J / if) dA/pp ratio as a function of rapidity and impact 
parameter [20] ■ The large cc total cross section also has implications for the J /if) yield 
if J/^'s arise from cc recombination in a QGP. Such a total cross section would suggest 
significant secondary J /if) production at RHIC, leading to enhancement rather than 
suppression in central collisions. There is no evidence for a strong regeneration effect in 
the PHENIX Au+Au data so far, see Thews' talk j2Tj . 

Shadowing, the modification of the parton densities in the nucleus with respect 
to the free nucleon, is parameterized as Ff(x, fi 2 ,b, z) = S l (A,x, p 2 ,b, z)ff(x, p 2 ) in 
Eq. ((H). We did not discuss the effect of shadowing on the charm pr distributions because 
the effect at midrapidity is small and, on the logarithmic scale of the pt distributions, 
negligible. The J /if) is another story due to the PHENIX muon capability at forward 
and backward rapidity. As shown in Pereira's talk 2\, although the PHENIX J /if) data 
are consistent with shadowing alone, the data are also consistent with nuclear shadowing 
plus a small absorption cross section of 1 — 3 mb, smaller than that currently obtained 
in SPS measurements (22] • We have calculated J / if) production in the color evaporation 
model (CEM) using the same mass and scale as in cc production but cutting off the 
invariant mass of the pair at Am 2 D . The calculations of the dA/pp ratios are done at LO 
to simplify the calculations. As shown in Fig- El the LO and NLO ratios are equivalent. 
We have now also implemented nucleon absorption in the calculation, showing the effect 
of several absorption and production mechanisms. 

To implement nuclear absorption on J /if) production in pA collisions, the pN 
production cross section is weighted by the survival probability, S ahs , so that |23] 

d 2 b / dzp A (b,z)S* hs (b,z) (4) 

J — oo 

where b is the impact parameter and z is the longitudinal production point. If S &hs = 1, 
a pA = Aa pN . For S abs ^ 1, a pA = A a a pN . We define S abs as 

S ahs (b, z) = exp j- J™ dz' PA (b, z')a &hs {z' - z)} . (5) 
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Figure 2. The J/ijj pAn/pp ratio at 200 GeV. We compare the NLO (solid histogram, 
MRST HO) and LO (solid curve, MRST LO) results using rn = /x/2 = 1.2 GeV with 
the EKS98 parameterization. 

The nucleoli absorption cross section, cr a b s , depends on where the state is produced 
and how far it travels through nuclear matter. The effective A dependence is obtained 
from Eqs. (J3J) and (JHJ) by integrating over z', z, and b. The contribution to the full 
A dependence in a(xp) from absorption alone is only constant if <r a b s is constant and 
independent of the production mechanism The observed J/ip yield includes feed 
down from XcJ and if)' decays, giving 

SfoQb, z) = 0.585^ dir (fo, z) + 0.3S£(b, z) + 0.12S$*{b, z) . (6) 

In color singlet production, the final state absorption cross section depends on the size 
of the cc pair as it traverses the nucleus, allowing absorption to be effective only while 
the cross section is growing toward its asymptotic size inside the target. On the other 
hand, if the cc is only produced as a color octet, hadronization will occur only after the 
pair has traversed the target except at very backward rapidity. We have considered a 
constant octet cross section, as well as one that reverts to a color singlet at backward 
rapidities. For singlets, Sfj^ Aix ^ ^ Sffi s but, with octets, one assumes that 

Sj% dir = ^Xcj = Sffi B - As can be seen in Fig. El the difference between the constant 
and growing octet assumptions is quite small at large \/S with only a small singlet 
effect at y < —2 and —5 at RHIC and the LHC respectively. Singlet absorption is 
also important only at similar rapidities and is otherwise not different from shadowing 
alone. Finally, we have also considered a combination of octet and singlet absorption 
in the context of the NRQCD model, see Ref. for more details. The combination 
of nonperturbative singlet and octet parameters changes the shape of the shadowing 
ratio slightly. The results are shown integrated over impact parameter for the EKS98 
shadowing parameterization since it gives good agreement with the trend of the PHENIX 
data. 

We will not discuss the spatial dependence of shadowing and absorption in any 
detail here. The spatial dependence of shadowing alone was discussed in Ref. j2D] • When 
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Figure 3. The J/tp dA/pp ratio with EKS98 at 200 GeV (left) and 6.2 TeV (right) as a 
function of rapidity for (a) constant octet, (b) growing octet, (c) singlet, all calculated 
in the CEM and (d) NRQCD. For (a)-(c), the curves are no absorption (solid), <7 a b s = 1 
(dashed), 3 (dot-dashed) and 5 mb (dotted). For (d), we show no absorption (solid), 
1 mb octet/ 1 mb singlet (dashed), 3 mb octet /3 mb singlet (dot-dashed), and 5 mb 
octet/3 mb singlet (dotted). 



absorption is included, the trend of the impact parameter dependence is in agreement 
with the PHENIX data at y > (the north muon arm) but is too weak to describe the 
strong dependence at y < (the south muon arm), see Pereira's talk 

3. Conclusions 

In summary, the RHIC d+Au data on open charm and J/ip are beginning to come 

into their own. While the QCD calculations agree well with the shape of the STAR px 

distributions, they underestimate the reported total cross section. In contrast, the J/ip 

cross section is in relatively good agreement with QCD predictions 2 and the agreement 

of the minimum bias data with calculations including shadowing and nucleon absorption 

is quite good. Work is ongoing to better understand the impact parameter dependence. 
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